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WING  DIVERGENCE  AND  ROLLING  POWER 
by 

LI.  T.  Niblett 

SUMMARY 

The  fundamental  static  aeroelastic  properties  of  a  swept  wing  which  has 
freedom  to  roll  are  studied.  It  is  shown,  using  a  simple  example,  that  the 
divergence  speed  of  a  wing  which  can  roll  steadily  is  high  or  non-existent 
and  consequently  the  rolling  power  of  any  aileron  will  vary  almost  linearly  with 
dynamic  pressure  and,  in  particular,  remains  finite  at  the  fixed-root  divergence 
speed.  Inboard  ailerons  have  higher  reversal  speeds  than  outboard  ailerons  when 
the  wing  is  swept  back  but  the  opposite  holds  true  when  the  wing  is  swept  forward. 
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I  INTRODUCTION 

During  an  investigation  of  the  effect  on  aileron  reversal  of  the  coupling 
of  flexure  and  torsion  of  swept  wings  through  stiffness'  it  was  realised  that 
the  properties  found  for  the  particular  case  of  no  stiffness  coupling  were  not 
all  that  familiar.  What  was  principally  in  question  was  the  significance  of 
divergence  to  rolling  power  at  practical  airspeeds.  Rolling  power  -  that  is 
rolling  velocity  divided  by  airspeed  for  unit  control  deflection  -  becomes 
infinite  at  a  divergence  speed  but  the  appropriate  divergence  speed  is  that  for 
the  case  in  which  the  rolling  moment  due  to  wing  distortion  is  reacted  by  the 
airforces  that  arise  from  its  rolling  velocity  rather  than  by  the  wing  root  being 

encastered  as  is  the  case  for  fixed-root  divergence.  Little  attempt  has  been 
2-4 

made  previously  to  estimate  the  speed  at  which  this  antisymmetric  divergence 
occurs  and  so  discover  its  practical  significance.  The  purpose  of  the  work 
reported  here  was  to  remedy  this  deficiency. 

The  previous  investigation'  did  not  include  rolling  power  as  a  subject  and 
that  omission  is  rectified  here  for  the  case  of  wings  with  no  stiffness  coupling. 
Arbitrary  modes  describe  the  structural  properties  of  the  model  as  they  are 
believed  to  be  the  most  efficient  and  instructive  fcrm  of  structural  representa¬ 
tion  in  cases  such  as  this.  First  the  general  equation  for  a  rolling  wing  is 
described.  From  this  equation  the  equations  for  fixed-root  and  antisymmetric 
divergence,  ie  divergence  when  the  wing  is  free  to  roll,  are  derived.  The 
rolling-power  equation  is  also  obtained.  The  rolling  power  can  be  evaluated  by 
dividing  one  function  of  speed  by  another.  The  zeros  of  the  dividend  give  the 
aileron  reversal  speeds  and  the  divisor  is  zero  at  antisymmetric  divergence 
speeds . 

The  equations  are  then  applied  to  an  example  wing.  This  wing  has  uniform 
spanwise  properties  and  is  assumed  to  distort  in  sets  of  modes 
1  -  cos(2m  -  I )  (tt/ 2) n  flexurally  and  sin(2m  -  1 )  ( tt/ 2) n  torsionally  (n  is  a 
non-dimensional  spanwise  coordinate).  The  aerodynamic  forces  are  based  on  strip 
theory  and  the  assumptions  that  the  aspect  ratio  is  large  and  the  sweep  small. 
Airforces  due  to  inboard  and  outboard  ailerons  am  represented  by  the  extremes 
of  discrete  forces  acting  at  arbitrary  chordwise  positions  at  the  root  and  the 
tip  of  the  wing.  Forces  at  the  root  neither  roll  nor  distort  the  wing  but  the 
ratios  between  the  different  actions  are  finite  and  from  them  can  be  deri'’°d  the 
mathematical  limits  of  the  benavicr  when  the  aileron  is  inboard  and  its  span 
tends  to  zero. 
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Initially  the  wing  is  taken  to  be  unswept  so  that  only  torsional  flexibil¬ 
ity  is  of  interest.  The  number  of  arbitrary  modes  included  in  the  calculation  is 
varied  to  find  the  minimum  needed  to  give  good  approximations  to  the  antisytaneiric 
divergence  and  aileron  reversal  speeds.  A  swept  wing  is  then  considered,  first 
with  only  flexural  flexibility  but  finally  with  both  flexural  and  torsional 
flexibility. 

2  THE  ROLLING  EQUATION 

2. 1  Description  of  equation 

The  equation  of  a  wing  '•oiling  with  constant  angular  velocity  can  be 
written 


_ 

— 

— 

—  —1 

C  +  E  u~2 
ss  ss 

b 

sr 

qs 

C 

SC 

c* 

rs 

b 

rr_ 

pi/v 

c 

rc  _ 

where  qg  is  the  vector  of  the  generalised  coordinates  of  a  set  of  arbitrary 
fixed-root  wing  modes,  Css  and  are  square  matrices  of  aerodynamic  and 

structural  stiffness  coefficients  and  is  a  scaled  dynamic  pressure,  c^.g  is 

a  row  matrix  of  aerodynamic  coefficients  derived  from  the  virtual  work  done  in 
roll  by  the  airforces  due  to  the  wing  distortions,  p  is  the  rolling  velocity  of 
the  wing,  i  is  a  typical  length  and  V  is  che  airspeed.  b^r  and  are  a 

column  matrix  and  a  scalar  derived  from  the  virtual  work  done  by  the  airforces 
due  to  the  rolling  velocity  of  the  wing  in  the  distortion  modes  and  in  roll  res¬ 
pectively,  Similarly  c  and  c  are  derived  from  the  virtual  work  done  bv 

sc  rc 

unit  displacement  of  the  ailerons.  For  a  rigid  wing  the  rolling  power  is  given 
by 


pr* 


V 


(2) 


2. 2  Fixed-root  divergence 

The  fixed-root  divergence  speeds  are  given  by  the  zeros  of  the  determinant 
|Css  +  E  o  |  .  A  wing,  prone  to  divergence,  probably  has  more  than  one  diverg¬ 
ence  speed  but  a  good  approximation  to  the  lowest  can  generally  be  obtained  using 
only  one  arbitrary  mode  each  of  flexure  and  torsion,  if  the  modes  are  chosen  with 
care.  The  fundamentals  of  the  modes  chosen  for  the  example  wing  whose  coef¬ 
ficients  are  evaluated  in  the  Appendix  are  t,  =  l  -  cos(ir/2)  i)  tut  flexure  and 
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0  =  sin(Tr/2)  n  for  torsion  and  Ic  1=0  if  there  is  no  more  than  one  flexure 

1  ss 

mode,  because  d£  /dq  =  (2m  -  1)(7t/2)  6  .  Hence,  if  the  flexural  axis  is  the 

m  m 

reference  axis  and  the  modes  are  not  coupled  by  structural  stiffness, 

-4  -4/ 


_2  —2 

Ic  +E  u  |  =  (c_,e  +c  ,e,;)u 

1  ss  ss  1  ff  tt  tt  ff' 


effettu  =  effettu 


t'-f)' 


(3) 


where  is  the  dynamic  pressure  at  fixed-root  divergence  and  is  given  by 


-2 

ude 


■  -  (—  * — V 

Veff  ett  / 


(4) 


Substituting  the  values  of  c„,  ctt>  eff  an^  erf  given  in  the  Appendix 
we  have 

(v)  °de  “  -((!)2a  ta"  A  (“)_1  +  «a(S)'')  (5) 

as  has  been  found  before* . 

From  this  we  see  that  aerodynamic  axes  upwind  of  flexural  axes  (?a  negative) 
and  sweep  forward  (A  negative)  lead  to  fixed-root  divergence.  The  first  of 
these  conditions  is  satisfied  subsonically  by  most  wings  but  if  they  are  swept 
back  they  do  not  suffer  fixed-root  divergence  if 

tan  A  >  -  (~)  A'*Ca|i  .  (6) 


2. 3  Antisymmetric  divergence 

Fixed-root  divergence  is  unlikely  to  be  a  physically-realisable  phenomenon 

outside  the  wind  tunnel.  Large  lift  forces  are  implied  in  the  symmetric  case 

and  these  would  lead  to  large  normal  accelerations  were  the  wing  part  of  a 

free-flying  aircraft.  It  has  been  suggested**  that  the  appropriate  case  is  one 

in  which  the  lift  force  is  reacted  by  the  aircraft  inertia  but  it  seems  far  more 

realistic  to  take  the  mode  of  free  divergence  to  be  one  in  which  the  incidence 

is  of  different  sign  inboard  and  outboard  so  that  any  incremental  upwards  lift  is 

6“  8 

opposed  by  an  equal  downwards  lift  .  This  gives  a  speed  about  twice  that  of 

fixed-root  divergence  but  this  does  not  mean  that  an  aircraft  can  be  flown  safely 

at  the  fixed-root  divergence  speed  for  an  aircraft  of  which  the  wing  is  a  part 

7  8 

becomes  difficult  to  trim  near  it  ’  . 
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-  c'  (c 
rs\ 


b  b_1c’ 
sr  rr  rs 


sU~2)  (C 


-  b  b-1c 
sr  rr  re/ 


Equation  (11)  shows  clearly  that  the  rolling  power  is  indeterminate  only  at 
the  antisymmetric  divergence  speeds  and  not  at  the  fixed  root  divergence  speeds 
and  that,  if  the  free  divergence  speeds  are  high,  rolling  power  varies  more-or- 
less  linearly  with  dynamic  pressure.  The  aileron  reversal  speed  can  be  found  by 
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equating  the  right-hand  side  of  either  equation  (10)  or  equation  (11)  to  zero  but 
a  more  convenient  and  informative  equation  is  obtained  by  rearranging 
equation  (I)  after  putting  p  zero 


r  v  "2 

C.  +  E  u 

ss  ss 

c 

sc 

1 

U1 

l_I _ 

0 

c ' 

_  rs 

c 

rc 

_  _ 

0 

_  _ 

which  can  be  reduced  to 

C  -c  c  *  c '  +  E  u  =  0 

|  ss  sc  rc  rs  ss  | 


(12) 


(13) 


an  eigenvalue  problem. 

3  UNSWEPT  WING 
3 . I  Antisymmetric  divergence 

When  the  wing  is  unswept  the  only  distortions  of  importance  are  those  in 

torsion.  The  modes  chosen  for  the  example  happen  to  be  exactly  those  of  fixed- 
.  -2 

root  divergence,  (C  +  E  v  )  is  a  diagonal  matrix  and  the  mth  element  on 
ss  ss 

the  diagonal  is 

i(>  -  (2m  -  1  )2»j/u2)  ,  (14) 

d 

2  . 

where  is  the  lowest  of  the  fixed-root  divergence  dynamic  pressures  as 

before . 

From  equation  (8),  with  the  values  of  the  coefficients  taken  from  the 
Appendix,  the  antisymmetric  divergence  pressures  are  given  by  the  solutions, 


If  the  series  is  truncated  to  only  one  term  the  solution  is  a  pressure  of 
2  7 

69.131)^  .  This  has  been  quoted  previously  as  an  approximation  to  the  free- 
divergence  pressure  but  more  than  the  fundamental  fixed-root  mode  is  needed  for 
an  adequate  description  of  the  wing's  distortion.  If  two  terms  of  the  series 
are  included  the  lower  of  the  two  solutions  is  8.1881;^  and  this  is  a  good 
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approximation  to  the  true  pressure  of  8.l83u^  .  The  wing  is  rolling  at  diverg¬ 
ence  and  Fig  I  shows  the  incidence  induced  by  the  roll  velocity  as  well  as  that 
due  to  distortion  of  the  wing.  The  net  incidence  resembles  sin(3ir/2)n  ,  the 

first  overtone  mode  of  fixed-root  divergence.  There  is  an  infinity  of  divergence 

2  2 

pressures  with  the  pressures  getting  closer  to  (2m  +  1)  as  the  order  m  gets 

higher. 

The  wing  is  uniform  in  the  present  case,  hence  0'  “  f  L  dn  ,  and  the 


angle  at  the  tip  in  any  of  the  divergence  modes  is  zero  for 
1  I  I  ) 


e 


3 0*j  9'dn  “  J  dri  f  L Cn j ) dn }  «  f 

0  On  6 


nL(n)dn  =  0 


0  On 

because  the  net  rolling  moment  is  zero. 

3 . 2  Rolli.ig  power 

The  extremes  of  possible  spanwise  positions  of  the  aileron  are  covered  by 
taking  as  examples  ailerons  whose  spans  tend  to  zero  (point  ailerons)  at  the  tip 
and  at  the  root.  When  the  values  of  coefficients  derived  in  the  Appendix  are 
substituted  in  the  right-hand  side  of  equation  ( ) 0) ,  its  zeros  are  given  by 


-  0 


m=l  ]l  -  (2m  -  I)2 


(16) 


for  the  tip  aileron  and 

'  -  « 


ra=  I 


(-1 )m* 1 ( 2m  -  1)  ' 
2\ 

2 

1  -  (2m  -  1)  -f 

u  , 


=  0 


(17) 


for  the  root  aileron. 

If  only  one  term  of  the  series  is  taken, equation  (16)  gives 

2 


(18) 


where  u  is  the  reversal  pressure. 
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If  another  term  is  included  a  quadratic  equation  is  obtained 
,4  -  x  2 


— -  -  (lO  -  7.38519 

\  U  „2 


+  1  -  0.900633  —  =  0 

a'  u  ^a 

r 


(19) 


The  solution  of  equation  (19)  can  be  put  in  the  form 


if 


18  ~  =  10  -  7.3851 

u  'a 

r 


9  ^  - /(8  -  7 . 385  1 9  -2. 


88202 


5:  10  -  7.38519 


H 


8  -  7.38519  ~ 

'a  \  ~*a 

x  (  1  -  I  .44101^ J\S  -  7.38519  ^ 


y 


^8  -  7.38519  ^ 


[8  -  7.38519  >  2.88202  -jr- 


(20) 


The  more  interesting  root  is  given  by 


2 

u 


1  -  0.820577 


0.010007 


(21) 


For  a  subsonic  wing  and  conventional  aileron  is  probably  somewhere 

between  -1  and  0  in  value  and  so,  in  the  present  context,  the  simplest  equation 
for  the  reversal  dynamic  pressure,  equation  (18),  gives  a  good  approximation  as 
it  agrees  with  the  two-term  formula,  equation  (21),  when  £  =  0  and  gives  a 

pressure  less  than  17  smaller  when  (C^/C^)  =  -•!  .  The  reversal  speed  is  below 
the  fixed-root  divergence  speed  if  the  centre  of  the  lift  due  to  aileron  rotation 
is  aft  of  the  flexural  axis. 


Substituting  the  simplest  approximations  to  the  divergence  and  reversal 
speeds  in  equation  (10),  modified  to  give  p  ,  the  rolling  power  of  the  flexible 
wing  as  a  fraction  of  that  of  the  equivalent  rigid  wing  (equation  (2)),  we  get 


The  ratio  is  almost  linear  with  dynamic  pressure  up  to  the  fixed-root 
divergence  speed  and  beyond. 


10 


The  reversal  speed  of  the  inboard  aileron  when  only  one  term  of  the  series 
is  included  is  given  by 


r 


and  when  two  tenns  are  included  an  approximate  lower  speed  is  given  by 


d 

~1 


1 .22608  +  0.047  I  59 

^a 


1.37934 


(24) 


The  reversal  pressure  given  by  the  one-term  approximation  is  only  about  32 
lower  than  that  given  by  the  two-term  approximation  when  fj./S  =  “I  and  the 
comments  made  on  the  tip-aileron  case  apply  equally  to  the  root  aileron.  Note 
that  the  reversal  speed  with  inboard  aileron  is  lower  than  that  with  an  outboard 
aileror  because  although  the  lift  of  an  inboard  aileron  is  a  poor  source  of 
wing  distortion  it  is  an  even  poorer  source  of  rolling  moment. 

4  SWEPT  WING 
4 . I  Antisymmetric  divergence 

It  has  been  shown  in  section  2.2  that  a  wing  with  fixed  root  and  only 
flexural  flexibility  diverges  at  some  airspeed  if  the  wing  is  swept  forward  but 
not  if  the  wing  is  swept  back.  However  if  the  wing  is  allowed  to  roll  it  does 
not  diverge  at  any  airspeed  whatever  the  sweep.  The  matrix  of  aerodynamic  stiff¬ 
ness  coefficients  for  the  flexural  modes  is  not  symmetric  because  the  airforces 
dep  ;nd  on  the  deflection  slopes  in  the  modes  rather  than  the  deflections  them¬ 
selves.  Hence  the  eigenvalues  associate,,  with  equation  (9)  need  not  be  real  and 
indeed,  if  enough  arbitrary  modes  are  caken,  all  are  purely  imaginary.  The 
eigenvalue  with  the  lowest  modulus  is 

u^a  =  (±40.75i)/lj  e1(A  tan  A)”'  (25) 


and  approximations  to  the  imaginary  constants  of  the  next  smallest  eigenvalues 
are  ±326i  and  ±11 OOi  .  The  rolling  power  of  a  wing  with  only  flexural  flexi¬ 
bility  is  therefore  never  indeterminate.  When  only  four  modes  are  taken  the 
dynamic  pressure  is  only  about  greater  in  modulus  than  the  true  divergence 
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pressure.  Purely-imaginary  eigenvalues  are  associated  with  rv.il  skew-symmetric 
matrices  and  this  is  consistent  with  the  number  of  degrees  of  freedom  needed  for 
a  good  approximation  being  about  twice  the  number  needed  when  the  eigenvalues  are 
real.  When  considering  the  differences  between  the  divergence  pressures  in  the 
free  and  encastered  wing  cases  it  should  be  remembered  that  in  the  free  case  the 
airforces  due  to  wing  distortion  are  only  reacted  by  the  airforces  due  to  rolling 
velocity  which  are  at  their  largest  at  the  wingtip  in  contrast  to  the  fixed-root 
case  in  which  the  distortion  airforces  are  only  reacted  at  the  root. 

<?  ,  the  variable  introduced  to  represent  the  r.  lative  amounts  of  flexural 
and  torsional  flexibility  of  the  wing,  is 


and  takes  some  of  the  aerodynamic  and  geometric  properties  into  account.  It  will 
be  assumed  that  is  negative  (tfj  is  positive)  which  is  not  unreasonable  fur 

a  conventional  wing  in  subsonic  flow,  and  hence  positive  values  of  ;■  represent 
sweptforward  wings.  Consideration  of  positive  as  well  as  negative  values  of  ^ 
might  increase  the  risk  of  confusion  without  adding  much  to  the  understanding  of 
what  happens  in  practical  cases. 

Fig  2  shows  the  variation  of  the  inverse  of  the  antisymmetric  divergence- 
pressure  with  flexibility  ratio.  When  the  flexibility  ratio  is  zero  the  diverg¬ 
ence  is  purely  torsional  and  it  can  be  seen  that  the  figure  is  limited  to  only 
the  lowest  four  divergence  roots.  Flexural  and  torsional  flexibility  seem  to  be 

independent  in  effect  and,  with  the  eigenvalues  identical  when  $  is  infinite, 

2  2 

whatever  its  sign,  the  curves  are  symmetric  about  the  u^/u  axis.  It  has 
already  been  noted  that  antisymmetric  divergence  has  an  insignificant  effect  on 
rolling  power  when  the  wing  is  unswept  and  it  appears  that  the  same  is  the  case 
when  the  wing  is  swept. 

A . 2  Rolling  power 

Rolling  peer  when  there  is  only  flexural  flexibility  is  easier  to  represent 
than  antisymmetric  divergence.  When  there  is  a  point  aileron  at  the  root,  c ^  is 
null  (see  Appendix)  and  tin  reversal  speeds  nr«  identical  uri  th  the  I  ixed-rtvt 
divergence  speeds  of  .ae-  wing  (equation  (  ’  )  Tin  two -node  approximation  leads 

to 


12 


I 


2  -2 
<  u 


-1 


0 


3 


i  +  8I< 


2  -2 
u 


(27) 


where  =  (tt/2)  "*EI(A  tan  A)  *  . 

From  equation  (27)  the  smaller 
cal  factor  of  which  can  be  compared 
(equation  (5))  and  -1.03968  for  the 
means  that  only  wings  that  are  swept 
at  the  root. 


reversal 
with  -I 
accurate 
forward 


pressure  is  -1.0406'  ,  the  numeri- 
for  the  one-term  approximation 
solution.  The  sign  of  the  factor 
suffer  reversal  if  the  aileron  is 


For  a  tip  aileron  the  two-mode 

-  ,  2  -2 
1  -  4/u  +  k  u 

3  -  4 /it 


approximation  leads  to 

-  (I  -  4/(3tt)  ) 

I  +  4/  (  3tt)  +  8K:2u~2 


0  . 


(28) 


From  equation  (28)  the  smaller  reversal  pressure  is  4.5039k  ,  the  numerical 
(actor  of  which  can  be  compared  with  3.66  for  the  one-term  approximation  and 
4.50958  for  the  accurate  solution.  The  sign  of  the  factor  means  that  only  wings 
that  are  swept  back  suffer  reversal  if  the  aileron  is  at  the  tip. 


Fig  3  shows  the  result  of  including  torsional  flexibility  as  well,  for 

three  particular  values  of  r  / f,  .  Mote  that  it  is  the  inverse  of  reversal 

b  3 

pressure  that  is  plotted  against  flexibility  ratio.  Again  low  reversal  speeds 
are  associated  with  tip  ailerons  when  the  wing  is  swept  back  and  root  ailerons 
when  the  wing  is  swept  forward. 


This  can  be  explained  using  simple  physical  considerations.  Imagine  an 
unswept  wing  with  an  aileron  which  has  been  rotated  trai 1 ing-edge  down.  The 
aileron  deflection  results  in  an  upward  lift  which  is  opposed  by  the  lift  due  to 
the  wing  twist,  which  is  a  consequence  of  the  aileron  lift  acting  aft  of  the 
flexural  axis,  and  by  the  lift  due  to  rolling  velocity.-  The  net  rolling  and  bend¬ 
ing  moments  must  be  zero  at  the  root  when  the  wing  reaches  its  equilibrium  roll¬ 
ing  velocity.  Thus  the  distribution  of  bending  moment  over  the  span  of  the  wing 
when  there  is  a  point  aileron  at  the  wing  tip  must  be  as  in  Fig  4a.  The  curve  of 
the  bending  moment  due  to  the  wing  airforces  must  be  concave  upwards  for  these 
airforces  are  one-signed  but  the  curve  of  the  bending  moment  due  to  the  aileron 
airforces  is  a  straight  line.  There  is  thus  a  sagging  bending  moment  over  the 
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wing  which  results  in  Li.  bending  tip  up.  Imagine  the  wing  to  be  swept  by  rota¬ 
tion.  If  it  is  swept  backward  the  lift  due  to  bending  is  opposite  in  sign  to  the 
aileron  lift  and  the  rolling  power  is  reduced,  although  there  is  some  alleviation 
of  this  because  the  wing  flexure  'blows  off',  but  if  the  wing  is  swept  forward, 
not  only  is  the  lift  due  to  bending  of  the  same  sign  as  the  aileron  lift,  but  it 
is  enhanced  by  the  lift  due  to  the  extra  wing  distortion  due  to  the  wing's  tend¬ 
ency  to  diverge  in  flexure.  When  the  aileron  is  inboard  the  net  loading  of  the 
wing  leads  to  a  hogging  moment  (Fig  4b)  which  has  quite  the  opposite  effect. 

The  dotted  lines  in  Fig  3  labelled  -R  indicate  the  negative  values  of 

2  2 

u./u  obtained  when  the  aileron  is  at  the  root  of  a  sweptback  wing  which  is 
dr 

flexible  in  flexure.  These  were  the  only  significant  negative  values  obtained. 

2 

If  consideration  is  restricted  to  cases  in  which  u,  is  positive  (£  negative) 

then  rolling  power  increases  with  speed,  for  this  particular  arrangement,  for  it 
2  2 

varies  nearly  as  I  -  u  /ur  .  Remembering  that,  on  a  rigid  wing,  an  inboard 
aileron  provides  less  rolling  power  than  an  outboard  aileron  because  its  moment 
arm  is  shorter,  there  appears  to  be  an  opportunity  of  optimising  the  rolling 
power  of  a  sweptback  wing  by  positioning  the  aileron  to  give  the  best  compromise 
between  basic  rolling  power  and  reversal  speed. 

It  can  be  seen  from  Fig  3  that  the  effects  of  flexural  and  torsional  flexi¬ 
bility  on  reversal  speed  are  almost  independent  of  each  other  if  the  flexibility 
ratio  is  anything  but  large.  Truncating  equation  (24)  and  using  the  exact  value 
for  the  flexure-only  case  we  get,  for  the  case  of  the  root  aileron  when  the  flexi¬ 
bility  ratio  is  low, 

Ud  5b 

— |  %  1  -  1.226  +  0. 96 24>  (29) 


remembering  that  $  is  positive  whet  the  wing  is  swept  forward;  and  similarly  in 
the  case  of  a  tip  aileron 


u 

u 


2 

d 

2 

r 


I  -  0.821  ~  -  0 . 2220 
a 


(30) 


That  part  of  Fig  3c  that  relates  to  a  root  aileron  duplicates  the  variation 
of  fixed-root  divergence  speed  with  scaled  flexibility  ratio  for,  with 
zero,  both  c^c  and  ctc  are  null. 
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Most  of  the  curves  were  found  from  four  flexural  and  four  torsional  modes 
as  data  but  checks  were  made  using  eight  of  each  type  of  mode.  The  values  calcu¬ 
lated  with  the  lesser  number  of  modes  are  good  approximations  to  the  lower  rever¬ 
sal  speeds,  the  more  important,  but  occasionally  are  not  as  good  for  the  upper 
speeds,  of  which  there  are  more  than  those  shown  in  the  Figure.  However  they  are 
always  accurate  enough  for  the  characteristics  of  the  curves  to  be  correct. 

h . 3  Ta i lerons 

Roll  control  is  obtained  on  some  aircraft  by  fitting  them  with  tailplanes 
whose  two  sides  can  be  rotated  differentially  in  pitch  (tailerons).  In  Lhis  case 
the  airforces  on  the  wing  are  not  influenced  directly  by  the  control  movements 
and  the  appropriate  rolling  equation  is  equation  (I)  with  c  null  and  c 
representing  the  rolling  moment  transmitted  through  the  fuselage  from  the  tail- 
eron.  The  fuselage  and  taileron  are  taken  to  be  rigid  in  themselves.  If  c 
is  made  null  in  equation  (11)  the  rolling  power  of  the  taileron  is  obtained  as 


b  pH 
rr 


l+c'(C  -b  b  *  c '  +  E  u  2)  *b  b' 

rs  ss  sr  rr  rs  ss  sr  rr 


(31) 


and  the  free  divergence  speed  is  again  the  speed  at  which  the  rolling  power  is 
indeterminate.  However  the  fixed-root  divergence  speed  is  also  significant  in 
this  case  for  if  Csc  is  put  null  in  the  aileron-reversal  equation,  equation  (13), 
it  becomes  the  same  as  the  fixed-root  divergence  equation,  it 


+  E  u 
ss  r 


=  0 


(32) 


Thus  the  variation  of  aileron-reversal  pressure  with  flexibility  ratio  is 
identical  to  that  for  a  root  aileron  whose  lift  acts  at  the  stiffness  axes  of  the 
wing  and  is  given  by  the  lines  marked  R  in  Fig  3c  and  rolling  power  can  increase 
with  speed  if  the  wing  is  flexurally  flexible  and  swept  back.  In  any  case  the 
variation  of  rolling  power  with  speed  will  be  near  linear  at  flight  speeds  as 
before . 

5  CONCLUDING  REMARKS 

A  study  of  the  rolling  equation  of  a  wing  has  shown  that  there  is  a  free 
antisymmetric  divergence  case  analogous  to  the  free  syimetric  divergence  case. 

When  symmetric  divergence  occurs  the  distortion  mode  is  such  that  the  net  lift  on 
the  wing  is  zero  and  when  antisymmetric  divergence  occurs  the  mode  is  such  that 
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the  net  rolling  moment  is  zero.  In  the  symmetric  case  the  wing  pitches  by  an 
amount  sufficient  to  equalise  the  upwards  and  downwards  lifts  on  i t  and  in  the 
antisymmetric  case  the  wing  rolls  at  a  velocity  sufficient  to  equalise  its  posi¬ 
tive  and  negative  rolling  moments.  In  both  cases  the  divergence  mode  is  overtone 
in  shape  and  hence  the  critical  speeds  are  much  higher  than  that  of  fixed-root 
divergence.  The  order  of  the  overtone  is  the  higher  in  the  antisymmetric  case 
and  so  the  divergence  speed  is  also  the  higher. 

Furthermore  only  wings  with  little  or  no  sweepback  have  real  free  anti¬ 
symmetric  divergence  speeds.  The  eigenvalues  corresponding  to  the  divergence 
dynamic  pressures  of  wings  with  only  flexural  flexibility  are  peculiar  in  that 
they  are  imaginary  numbers.  This  means  that  such  wings  do  not  diverge  and  the 
eigenvalues  are  seen  to  be  so  large  that,  at  practical  airspeeds,  rolling  power 
is  little  affected  by  possible  divergence-type  distortion  of  the  wing.  When 
torsional  flexibility  is  included  the  dominant  eigenvalues  change  first  to  con¬ 
jugate  complex  numbers  and  so  the  wing  still  does  not  diverge.  The  behaviour  of 
the  wing  is  identical  for  equal  amounts  of  sweepback  and  sweepforward  (Fig  2) 
and  this  is  consistent  with  the  aerodynamic  coefficients  for  lift  due  to  torsion 
having  the  same  values  for  equal  amounts  of  sweepback  and  sweepforward  and  the 
coefficients  for  lift  due  to  flexure  being  equal  but  opposite  in  sign  in  the 
same  circumstances. 

It  is  the  free  antisymmetric  divergence  speed  and  not  the  fixed-root  speed 
that  could  be  of  significance  in  the  calculation  of  rolling  power  but  as  it  is 
high  or  non-existent  it  is  of  little  consequence  in  practice.  If  the  airforces 
are  linear,  rolling  power  is  very  close  to  being  a  linear  function  of  dynamic 
pressure  being  zero  at  the  aileron  reversal  pressure.  This  reversal  pressure 
depends  on  the  spanwise  position  of  the  aileron  as  well  as  the  sign  of  the  sweep 
of  the  wing.  It  is  well  known  that  a  wingtip  aileron  on  a  swept  back  wing  has  a 
low  reversal  pressure  but  the  reversal  pressure  for  an  inboard  aileron  on  the 
same  wing  can  be  negative  which  means  that  rolling  power  increases  with  dynamic 
pressure  although  its  absolute  value  might  be  small  owing  to  the  shorter  moment 
arm. 

The  calculations  demonstrate  the  economy  that  can  be  achieved  when  well 
chosen  but  simple  arbitrary  modes  are  used  for  the  description  of  structural 
distortion  in  studies  such  as  this.  In  nearly  all  instances  two  such  modes  of 
each  type  have  been  shown  to  be  enough  to  give  results  sufficiently  accurate  for 
practical  purposes. 
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EVALUATION  OF  THE  COEFFICIENTS  FOR  AN  EXAMPLE  WING 
A. I  Example  wing 

The  example  wing  is  of  uniform  chord  c  perpendicular  to  the  stiffness 

axis  which  is  of  length  £  and  sweepback  A  .  Air,  of  density  p  ,  is  flowing 

past  the  wing  at  speed  V  .  Let  the  downward  heave  of  the  stiffness  axis  in  the 

flexural  mode  (qr  )  be  C  c  ,  where  4  =  I  -  ccs(m  -  !)nri  and  r,  is  zero  at 
fm  m  m 

Che  wing  root  and  unity  at  the  tip,  and  the  nose-up  incidence  (sections  perpen¬ 
dicular  to  the  stiffness  axis)  in  the  torsion  mode  (q  )  be  9  ,  where 

tin  m 

9  »  sin(m  -  l)nq  .  The  mode  appropriate  to  roll  (q  )  is  rotation  about  the  axis 

m  r  r  r  t  j- 

in  the  direction  of  the  line  of  flight  which  passes  through  the  wing  root.  Let 
the  generalised  coordinate  associated  with  aileron  angle,  taken  as  constant  along 
the  aileron  span  and  positive  trailing  edge  down,  be  q 

A. 2  Aerodynamic  coefficients 

The  work  done  in  a  small  displacement  of  the  wing  is 


w  -  -  Y  j 

L(c4i  qf  +  £n 

cos  A  dq^)  + 

/ 

M(86qt 

-  sin  A 

sqr> 

wing 

wing 

=  pc 2£V2  cos  A  a  jfiq^ 

fqt  £q;] 

"cff 

Cft 

b  ' 

fr 

qf 

,  (A- 1 ) 

ctf 

Ctt 

b 

tr 

qt 

_crf 

crt 

b 

rr_ 

(p£/V)q 

where,  following 

9 

Yates  , 

L  (upward)  = 

2  2 

-  pc£V  cos 

A  a(YC'qf  + 

9qt  + 

(pl/V) 

cos  l 

\  nq  , 

t+  Vc} 

,  (A-2) 

2  2  2 

M  (nose  up)  =  pc  £V  cos  A  a  {£  (yC'qf  +  9q,  +  (p£/V)cos  A  nq  )  +  £,  £0q  ),(A-3) 

a  I  C  V  D  D  C 


where  a  is  the  dC./da  of  the  wing,  i.  multiplied  by  a  is  dC,  /dq  and 
L  8  L  c 

and  are  the  chordwise  distances  of  the  wing  and  aileron  lift 

respectively,  aft  of  the  stiffness  axis  and  scaled  by  the  chord.  4'  =  d^/dn  and 
y  =  A_1  tan  A  where  A  =  i/c  and  it  has  been  assumed  that  the  value  of  y 
is  low,  ie  the  aspect  ratio  is  high  and  the  sweep  is  low,  and  so  terms  which  have 
as  a  factor  are  negligible  compared  with  unitv. 
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The  ailerons  considered  can  be  thought  of  as  point  ailerons  because  their 
spanwise  dimension  is  taken  to  be  almost  zero.  They  are  taken  as  located  at 
cither  the  root  or  the  tip  of  the  wing.  Only  ratios  of  the  coefficients  are 
required. 


A  cos  A  c^  cCrc  ~  -  cos(m  -  $)7Tn)/n  -*-0  ,  n  -*•  0 


1  ,  n  -  1 


,  (A- 5) 


A  cosA  ctn  c/crc  =  sin(m  -  i)irn/n  -*■  (2m-  1 )  (tt/ 2) ,  n  -*•  0 


(  1 )  »  n  —  1  • 


A.  3  Structural  coefficients 

Twice  the  strain  energy  of  the  wing  in  an  arbitrary  distortion  is 
I  1 

2V  =  (A2Z)~'eiJ  C"2q2dn  +  i_lGJ  j  6,2q2dn  »  pc2*V2a  [q*. q^J  |Eff  Efc  J  (qf, qt }  (say), 


•  (A-6) 


where  is  a  reference  speed. 


pc2J.V2ae,  ,  »  (A2*)  1  El  [  c"?"dn 

0  fm,fn  J  smsn 


1 

(2m-  l)2(2n-  I  )  2(tj/2)^(A2£)  1  El  J  cos(m  -  {)irn  cos(n  -  })irndri 

0 


=  J(2m  -  1  )t*  (it/ 2 )  ^(A2i)  *EI  ,  m  »  n  ;  =0,  m  i  n 


(A- 7) 


Similarly 


2.„2 

pc  £V„aet 
p  0  tm,tn 


£-lGJ 


I 

[  e’e 

I  mn 


dq  =  j(2m  -  1)2(u/2)2£  *GJ  ,  m  =  n;-0,  mj*n 
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E  is  a  diagonal  matrix, 

ss 

pff’  Etcj  =  -  l)W2)4EI/A2,  ((2m  -  1)tt/2)2GjJ  ,  (A-8) 

2  2  2  ^ 

where  El  =  pc  Z  V^aEI  ,  etc. 

The  scaled  airspeed  is  obtained  from  the  relationship  between  the  different 
scaling  factors  used  for  the  aerodynamic  and  structural  stiffness  coefficients 
(equations  (A-3)  and  (A-7))  as 


V  cos  A 


and  a  typical  form  for  the  divergence  equation  is 


_  o 

|c  +  Eu  |  = 


(A-9) 


0 


(A- 10) 
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